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Fields in Electrically Short Ground Systems: 
An Experimental Study 

A. N. Smith 1 and T. E. Devaney 1 

(April 20, 1959) 

An experimental study of magnetic field distribution in a simplified radial ground 
system on poorly conducting soil under an electrically short, toploaded monopole is described. 
It is shown that the distribution is that expected from the theory of J. K. Wait in those 
portions of the radial system satisfying the assumptions of the theory, and that the theory 
may still be successfully applied for //-field power loss computations even when this is not 
fully the case. The particular model system studied exhibits a condition suggesting damped 
standing waves on the radials in the area where the radial spacing exceeds that required 
by the theory. 



1. Introduction 

For an electrical^ small monopole, the ground- 
system power losses very often are a controlling 
factor in radiation efficiency. The decision to con- 
struct at Cutler, Maine, a new high-power, very-low- 
frequency station with specified efficiency greatly 
exceeding that of any such existing station has thus 
resulted in a large amount of recent interest in 
ground-system design. 

The final ground-system configuration for Cutler 
came basically from computations based on Wait's 
[1, 2, 3] equivalent surface impedance formula for 
radial current (//-field) power loss for a plane, radial 
grid of wires laid with perfect contact on an earth- 
air interface; and from design modifications brought 
about according to ground charging current (jE-field) 
power loss computations based on formulas due to 
Abbott [4], Wheeler [5], and Wait [6]. The sound- 
ness of the theoretical design was partially and in- 
directly confirmed at the l.S. Navy Electronics 
Laboratory, San Diego, Calif., by comparing losses 
derived from impedance measurements [7 J with 
those computed from theory for a model toploaded 
vlf monopole over a radial-wire system and also over 
a perfectly conducting ground [8]. A somewhat 
more direct test of theory- was made by W. G. Hutton 
and others of Smith Electronics, Inc., Cleveland, 
Ohio [9, 10], who measured the field distribution (mag- 
nitude only) in a full-scale, fan-shaped ground-S3^stem 
sector constructed at the Cutler site. An analysis 
of his data by Wait [11] showed that the measured 
field distribution did indeed generally agree with 
that expected from a division of ground system 
current according to the surface impedance ratio 
derived from the theory. 

However, there existed sufficient uncertainty about 
some conditions in all these experiments that a re- 
fined field-distribution study seemed desirable. 
Accordingly, a phase- and magnitude-measurement 
program was carried out on the vlf model previously 
used for impedance measurements. The dual ob- 
jective was to test the validity of the theory and to 
investigats more adequately the direction and degree 
of departures therefrom those existing in the model. 

'U.S. Navy Electronics Laboratory, San Diego, Calif. 



2. Theoretical Discussion 

The working formula for ground loss li g in terms 
of the azimuth al magnetic field Il^prf) at the surface 
of the ground (0=0) referred to the antenna feed 
point where current I Q flows is 



iL=Real part of 4 f Z[n 4> { P fl)] 2 2>Kp<Ip. 



(1) 



This equation is derived [2, 3] from the incremental 
Poynting vector flux due to the tangential component 
(E p ) of electric field arising from the finite conduc- 
tivity of the soil. In the derivation use is made of 
an approximate boundary condition through which 
the flux is expressed in terms finally of 11$, the mag- 
netic field of the monopole over an infinitely con- 
ducting ground. For convenience the integration is 
performed over two regions, ri<p<r, and p>r, r 
being the length of the wire radials, since, in the 
area occupied by the ground system, the surface 
impedance Z is the combination of wire impedance 
Z w and soil impedance Z s and in the remainder it 
is simply Z s . r t is the inner extremity of the radials 
The five major assumptions made in deriving equa- 
tion (1) are stated below. 

(1) The approximate boundary condition K p =ZH <P 
is used, requiring that pH change slowly in a distance 
d=\y~ 1 \ where y=^[ifxo)(a-{-io:ke )] y2 is the propaga- 
tion constant of the soil. 

(2) The combined surface impedance varies slowly 
in Jt" 1 !. 

(3) The radial wire screen of N radials in any 
small region a distance p from the antenna base 
equals that of a grid of parallel wires of the same 
mean spacing s=2wp/N, i.e., Z w =(iri s/\) In (s/2r)c) 
where rj =S77 ohms, and c is the wire cross section 
radius. This requires the condition s<C|Y e _1 | where 



f yl+iy 2 



is the effective propagation constant 



of the composite medium at the interface (y = iP)- 

(4) To replace the actual magnetic field H^ by H^ 

at the surface of the ground it is necessary that 
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It -1 !^!//^ i.e., the soil skin depth be much less than 
a free-space electrical radian. 

(5) The surface impedances of screen and soil 
must combine as simple parallel impedances, i.e., 
mutual coupling and potential differences between 
wires and soil are zero, so that 



Zj — Zj s Zj w j \Zj s -\- Zj w ) 



(2) 



within the ground system. 

These assumptions must be satisfied in setting up 
any direct experimental test of equation (1). 

Since this paper is largely concerned with the field 
structure due to division of current in the ground 
system according to assumption (5), its implications 
are now developed. The total ring current crossing 
a cylindrical boundary of radius p drawn about the 
monopole base as a center is I t = T w -\-I s =2TpH (l) . 
Normalizing with respect to base current (A =^// o ) 
and removing the factor 2irp, assumption (5) leads 
to 



h^ — h^Z s l\L w -\-Z^) -\-h<f>Z w j \Z w -\- Z s ) . 



(3) 



The first term on the right of (3) corresponds to 
wire current I w . Letting I w /I = Ni w , so that i w is 
the normalized current in one radial wire, the I w 
term is seen to be exactly Ni w /2iirp. Hence the mag- 
netic field h w at a field point d meters from the wire 
center is 



h w =h 4 ,pZ s /Nd(Z w +Z s ) . 



(4) 



The second term on the right of (3) corresponds to 
current in the soil I s , and in fact is i s /2irpj where 
i s =I s /I . If the wire centers lie exactly in the plane 
of the earth-air interface, then the magnetic field 
tangent to it, h s , can be simply evaluated by applying 
Ampere's law to the path illustrated in figure 1 , where 
n, the number of skin depths 8 from the surface to the 
portion def of the path, is so large that H^ along it is 
negligibly small. At the interface there is no tan- 
gential contribution from the wires themselves, and 
along the vertical portion the contributions cancel. 
Hence 



2T P h 8 =Ni w /2 +i s . 
Subtracting this from (3), it is seen that 



K ~2 \} + Z,+zJ 



(5) 



(6) 



That is, in the plane of the wire centers the total 
tangential field is the incident field of the downlead 
plus a contribution due to the induced current flowing 
in the soil. The wires carry only a portion of the 
image current. In the special case of a wire counter- 
poise and no soil, the second term vanishes, and all 
the image current flows in the wires. In the case of 
infinitely conducting ground, on the other hand, the 
second term also equals hJ2 (perfect reflection), the 




Figure 1. Integration path for evaluating h s . 



distinction between h w and h s is lost, and the total 
field is just hf. 

From the above considerations one sees that to 
determine the validity of the theory of 77-field power 
losses, it is sufficient to observe that the fields h w and 
h s satisfy equations (4) and (6) when assumptions 
(1), (2), "(3) and (4) are met. 

3. Experimental Procedure 

In the Cutler design [11] and in the NEL models 
used for the impedance measurements, soil constants 
and ground system parameters were such that as- 
sumptions (1), (2), and (3) were certainly met except 
at points either very close to the downlead or within 
a soil skin depth of the radial ends. For the purpose 
of the present experiment, condition (4) was regarded 
as an object of interest, while condition (5) was the 
main point at issue. 

Preliminary computations indicated that there 
would be sufficient individual wire current to yield 
usable signal for the instrumentation employed only 
if relatively few radials made up the ground system. 
This condition was required also to insure that the 
signal from a field point over a wire would be large 
compared to that from the nearby soil current, and 
moreover to insure that the wire field would not con- 
tain significant contribution from adjacent wires. 
To a radius of 100 ft the soil surface to a depth of 
several inches was stripped of all old radials, vegeta- 
tion, and inhomogeneities, and was regraded. Eight 
No. 30 bare copper wires 40 ft long were laid, 
stretched, and covered with a thin layer of mud to 
make contact with the interface as ideal as possible. 
The inner termination was a 1-ft-radius brass plate 
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grounded to six 18-in. copperweld rods; the outer 
terminations were % in. by 6 in. galvanized boat 
nails. Later, another sel of eight radials was run in 
between the originals. Finally all were extended to 
80 ft. The antenna was otherwise left as it had been; 
it is a somewhat conical, 36-radial topload supported 
at its periphery by a hexagonal spreader stretched 
over six 8-ft insulated galvanized iron pipe towers 
symmetrically spaced on a 13-ft-radius circle; the 
central support is the 12%-it downlead. At 1.5 Mc, 
the only frequency employed for both phase and 
magnitude measurements (an extensive series of 
magnitude distribution studies was made at 3.0 as 
well as 1.5), this structure on an ideally conducting 
ground has a base impedance of (0.40-J290) ohms, 
and has an effective height of 7% ft. 

Figure 2 shows the measurement scheme. It 
comprised two 1-in. -radius, shielded, 20-turn pickup 
coils Ci and C 2 connected by their coax lines and 
trimmers T x and T 2 through similar d-c amplifiers 
Ai and A 2 to the deflection system of an oscilloscope. 
Coil Ci was kept at the reference point 1 ft from the 
downlead over a radial wire, while C 2 was the ex- 
ploring coil. Both coils were oriented for maximum 
coupling, and the shields of both were insulated 
from direct contact with soil or wires. An RC phase 
shifter could be introduced into the circuit of one of 
the coils in order to sense the absolute direction of an 
observed phase shift. Base current was monitored 
by a calibrated milliammeter. The sensitivity of the 
field-detection system was obtained by comparison 
with readings on a field-strength meter (f-s meter- 
No. 1) set up beyond the radials, and also by meas- 
urement of the induction field of a known, long, 
straight current filament. Harmonic content of the 
antenna current was checked as a part of the cali- 
bration procedure. 

Since the distance d from wire to field point was 
of the same order as the loop diameter (d=1.94 in.) 
the field cannot be considered uniform over the coil 
area. The average field h w actually measured is re- 
lated to the value at the field point by a factor 




h w =Fh w where for coplanar radial and coil turns 



Figure 2. Scheme for measuring field magnitudes and phas 



F= 



b\'l-b 2 /d 2 



'2d 2 (l—y/l-b 2 /d 2 ) 



(7) 



For the coils used this amounted to 0.80. Calibra- 
tion against f-s meter No. 1 gives h w =KS w /I , where 
S w is the deflection developed on the oscilloscope, 
and K is the calibration figure. Stability of K 
turned out to be to within a percent over an ex- 
tended period. 

For measuring the field between the wires, small 
trenches were dug and the coil centers alined in the 
plane of the wires within one-fourth in. These 
measurements were done largely along the bisectors 
of the angle between adjacent wires, but explora- 
tion was also carried out azimuthally for several 
fixed values of p to test for constancy of field be- 
tween wires. 

The test of equivalence of 11$ and 11% was a 
comparison between the near-zone field measured 
for the antenna over the brass center plate only and 
the values computed from Wait's formula [12] 



^=%r [Gqv)' /2 +^ ta,1_1 H +z? 



(8) 



where D is the contribution to 77 from the topload 
disk, and h is the height of the disk above the ground 
(also taken to be the effective height). This was 
done using an untuned probe with a second field- 
strength meter calibrated against f-s meter No. 1. 
Later a partial confirmation was had from a top- 
loaded monopole set up over a small ground system 
half of which was hardware-cloth mesh, the other 
half consisting of 25 radials. 

Finally, it was necessary to know something about 
the electrical characteristics of the soil. To deter- 
mine these, the properties of short, known lengths 
of buried open-wire transmission lines of several 
spacings were measured at the frequencies of interest. 
From this information the electrical constants of 
the medium were derived; the best values for 1.5 
Mc appeared to be conductivity a=0.010 mho/m, 
dielectric constant k = 50. 

4. Discussion of Results 

Figure 3 shows the results for the test of assump- 
tion (4). In general, the agreement between the 
measured h<f> and theoretical h% is good. The de- 
parture within the periphery of the topload is not 
unexpected inasmuch as (a) the topload was some- 
what conical rather than flat as assumed in com- 
puting h% , and (b) the soil constants imply a skin 
depth of about 12 ft at 1.5 Mc, so that at no point 
under the topload can the currents be regarded as 
flowing in a plane, radial sheet without vertical 
divergence, which has the effect of decreasing the 
apparent field. The increase in h$ over the theoreti- 
cal curve beyond the topload seems to be mostly 
due to an 8 percent calibration error discovered in 
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Figure 3. Comparison of measured and theoretical h0 for 
homogenous, poorly conducting ground. 



f-s meter No. 1 after the conclusion of the experi- 
mental program. Thus for soils no poorer, 
relatively speaking, than that in this experiment, 
confirmation seems to be had of the propriety of 
replacing h^ by h$ in ff-field power loss computa- 
tions for the region beyond the ground system. 

Analysis of the radial ground system experiments 
proceeds from the observation that there seems to 
be only qualitative agreement between theory and 
experiment, and a significant difference between 
phenomena in the 40-ft system and those in the 80-ft 
(see figs. 4, 5, and 6). It is first in order to deter- 
mine the degree to which these somewhat sketchy 
ground systems satisfy the assumptions of the theory. 

Since pH changes at most by a factor of two in 
the ring 8<p<28 {8 being the skin depth) and less 
rapidly thereafter, the first assumption appears to 
be met fairly well except for the area under the top- 
load. The second assumption also appears to be 
satisfied for the same range of p for the 16-radial 
ground system, but not for the eight. The third as- 
sumption was not satisfied by the 80-ft 16-radial 
system for p> 40 ft; only poorly in 20 ft <p< 40 ft, 
for both 40- and 80-ft systems; fairly well in p< 
20 ft. In the eight-radial system it was not satisfied. 
Assumption (4) seemed to be well satisfied not only 
from the fact that at 1.5 Mc \y\ = 10/3, but also from 
the experimental confirmation discussed above. 
Finally, assumption (5) seemed to be only approxi- 
mately satisfied. Azimuthal sweeps at constant p 
over a quadrant of the ground system showed con- 
stancy in \h s \ within a few percent for all locations 
except those within a horizontal distance of 0.1 8 to 



£ 

■e- 



-4-H-H 


1 I 1 1 - 








~ 






HOMOGENEOUS 






h w 


1 MEASURED 






"*^^ys 


2 


- V 








_ 


V \ 


1 




- 




\ 








_ 


\ \ 4 


Z 


- 


\ ^ ft. 


- 


1 1 1 1 


I I i I ^ 



50 100 

/°. ft. 

Figure 4. Measured field magnitudes over and between wires 
of 40-ft 16-radial system. 
Short-dashed curve is \h w \ — \h,\, representing minimum possible h w . 




Figure 5. Comparison of measured and theoretical magnitudes 
over and between wires of 80-ft 16-radial system. 

Dash-dot curve is empirical itc=Aey x +Be—v*. 
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Figure G. Comparison of measured and theoretical phases for 
40- and 80-ft 16-radial ground systems. 
Positive phase represents lead relative to downlead current. 

the wires, where it showed a decided decrease. This 
one should expect from the interaction of fields of 
the parallel-flowing wire and soil currents, but the 
effect is quite localized. 

With reference to the above, then, for p< 20 ft 
measured \h s \ should be expected to fall below the 
computed values. Figures 4 and 5 show this to be 
true for both 40- and 80-ft systems. However, for 
p> 20 ft, the behavior of \h s \ is distinctly different, 
the 40-ft system showing reduction in \h s \ at all 
points, while the 80-ft system shows an enhancement. 
In the 40-ft system, the wire field in the first 10 ft 
is very nearly that predicted by current division 
according to parallel impedances, while for the 80-ft 
system agreement is only fair. In both systems, 
beyond 20 ft an increasing departure between meas- 
ured and computed \h w \ is evident, the nature of 
which appears different for the two situations. 

It is tempting to say simply that, since beyond 
20 ft the spacing s between the wires violates the 
condition s^^lTe" 1 !? they are acting no longer as a 
grid but as a collection of N independent, parallel, 
current filaments, each of whose internal impedance 
Zi=Ri-\-jXi is given by the usual formulas [13]. 
Such a radial wire mesh would have at a distance p 
from the common center an equivalent surface im- 
pedance of Z w =2irpZi/N. At 40 ft, the sixteen No. 
30 - copper - radial system would by this have 
Z w =2. 89 /38.1° . Since the soil impedance in parallel 
is Z,=33.1 /32.7° for (7=0.010 and k = 50, eq (4) 
would yield h w = 0.0475 m" 1 for this distance; for 
p=60 ft, h w = 0.022 m" 1 . Experimentally, one finds 
in the 40-ft system A w = 0.055 at 39 ft; in the 80-ft 



system, 0.052 and 0.0183 at 40 and 60 ft, respec- 
tively. Such agreement seems good until one ex- 
amines the relation between observed and predicted 
phases, particularly in the light of the contrasting 
behavior of \h s \, for the two radial lengths, in the 
region p> 25 ft. Thus, it would seem that the 80-ft 
s\^stem acts not merely as a grid of independent 
wires in parallel with the soil in the outer two thirds. 
The nature of the departure of \h w \ from the theo- 
retical values is somewhat suggestive of standing 
waves. Assuming a transmission-line mode of be- 
havior, one represents the current i w in the wire at a 
distance x from some reference point by [14]. 

iu,=Ae y e x +Be- y ^ (9) 

where y e is the effective propagation constant of the 
medium surrounding the wire, and the ratio B/A is 
the reflection coefficient. Regarding the outer ter- 
mination as the origin and looking back toward the 
downlead as the generator, the empirical fit with 
the experimental curve illustrated by the dash-dot 
curve of figure 5 is obtained for 7 = 0.10+^0.197, 
^4—1, B= — %, i.e., a lossy line terminated in an im- 
pedance higher than its characteristic impedance. 
Confirmation is had from the figures for a and k 
derived from y e 2 = ( 7o 2 +7 2 )/2. These are 7 = 0.0064 
mho/m, and k = 56A; also \=32 m, and ?;=0.16 c 
where c is the velocity of light in a vacuum. The 
buried -transmission-line measurement gave respec- 
tively 0.010, 50, 30, and 0.15, which is a fairly reason- 
able agreement. Since in actuality the wires were a 
little below the interface, a and k probably were be- 
tween the above and the values 0.0033 and 29, re- 
spectively, obtained from y=y e (completely buried 
wire) . 

The observed phases seem to fit generally with the 
above picture. The change should be 180° in 
X/2=16 m. The measured wire field phase referred 
to the downlead does not change* by 180° in this 
distance, but advances more slowlv at the rate of 
1.90°/ft instead of 3.40°/ft. However, if at each 
point one considers the total phase change of wire 
field with reference to that at the same radius, 
the combined change is 1.86°+1.29° = 3.15°/ft, out 
to about 50 ft where the total difference is about 
130°. The return toward zero for the phase of 
the field between wires for p>50 ft may be due to 
masking of the standing-wave circulating-current 
field in the soil by the incident downlead field. In 
the case of the 40-ft system, some evidence of the 
same thing shows up in the phase progression where 
in the soil it is 2.5°/ft, in the wire 1.3°/ft, total 3.8°/ft, 
but the magnitude curves, figure 4, show no recog- 
nizable standing wave characteristics. This may be 
because the region of independently acting radial 
length is considerably under half a wavelength (for 
the wire in the interface). 

5. Conclusions 

In both the 16-radial systems as well as in the 8- 
radial, it appears that a substantial potential differ- 
ence existed between the wires and the soil. Hence, 
in no case studied experimentally can the grid and 
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soil be considered as surface impedances in simple 
parallel combination. However, in those regions 
where the ground-system parameters satisfy the 
remaining four assumptions made in the theory, the 
observed fields seem to agree with the predicted ones 
fairly well. In view of the similarity of soil condi- 
tions for the present experiment and the previous 
impedance studies on systems of many radials, it is 
likely that conclusions drawn here hold for those 
earlier studies as well. 

To compute the iJ-field losses in a radial ground 
system laid on poorly conducting soil, one properly 
should consider four regions separately. The inner 
region is that in which the fields vary more rapidly 
than the limit specified in the theory; the second (in 
a proper design, by far the largest) is that in which 
at least the first four assumptions are satisfied; the 
third, that in which the wire spacing is such that 
the radials act independently; the fourth, the termi- 
nal region within a soil skin depth of the radial ends. 
In all four regions, copper losses would be correctly 
computed only from an observed current distribu- 
tion, but it seems likely that an adequate approxima- 
tion would result from using the distribution com- 
puted from the parallel impedance assumption; in 
any reasonable case these losses would be small. The 
ground losses for the inner region would have to be 
obtained from a volume current density distribution; 
they could be (usually are) eliminated by using an 
effectively solid sheet of highly conducting material 
in place of the radials. In the second region, where 
all assumptions but the fifth hold, the losses should 
be computable from Wait's theory even with very 
poor contact between wires and soil. The greater 
the number N, the better the results should be for 
losses computed in this second region. In the third 
and fourth regions, one has recourse to numerical 
integration of J*Z s H s 2 dA where H s (p, 0) would have 
to be an experimentally determined quantity. For 
very poor soil, poor contact between radials and soil, 
and mismatched radial termination, sufficient stand- 
ing-wave circulating current may exist to add a sig- 
nificant extra iJ-field loss component. Under these 
conditions it would seem worthwhile to direct atten- 
tion toward providing a termination in the character- 
istic impedance of the radial regarded as a transmis- 
sion line. 

It must be remembered that the above-described 
conditions are also those for which considerable 
charging current power loss can appear in the contact 
impedance between radials and soil, as well as in the 
displacement flux terminating in the soil between 
radials, and in fact such "E-fLeld" losses can be the 
dominating term in the total ground system loss 
picture, particularly where the loss tangent of the 
soil is of the order of unity. For better soils, these 
losses become of less importance, along with the 
standing-wave phenomenon, until for <r/uke o y>^>10 
they likely are of little significance. 
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